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The variational cluster approximation is used to study the frustrated Hubbard model at half filling 
defined on the two-dimensional square lattice with anisotropic next-nearest-neighbor hopping parameters. 
We calculate the ground-state phase diagrams of the model in a wide parameter space for a variety of 
lattice geometries, including square, crossed-square, and triangular lattices. We examine the Mott metal- 
insulator transition and show that, in the Mott insulating phase, magnetic phases with Neel, collinear, and 
spiral orders appear in relevant parameter regions, and in an intermediate region between these phases, 
a nonmagnetic insulating phase caused by the quantum fluctuations in the geometrically frustrated spin 
degrees of freedom emerges. 


1. Introduction 

The effect of geometrical frustration in strongly cor¬ 
related electron systems has been one of the major is¬ 
sues of condensed matter physics. In particular, a spin- 
liquid state caused by the frustration has been inter¬ 
preted as an exotic state of matter, where the magnetic 
long-range order is destroyed, yielding a quantum para¬ 
magnetic (or nonmagnetic) state at zero temperature^^ 
or even exotic mechanisms of high-temperature super¬ 
conductivity.^^ The Hubbard, Heisenberg, and related 
models defined on two-dimensional square and triangular 
lattices with geometrical frustration have been studied in 
this respect to find novel quantum disordered states by 
a variety of theoretical methods. 

In the square-lattice cases, the J 1 -J 2 Heisenberg 
model with the nearest-neighbor (Ji) and next-nearest- 
neighbor (J 2 ) exchange interactions have been stud¬ 
ied for more than two decades.At J 2 = 0, 
where the frustration is absent, the model is known 
to have the Neel-type antiferromagnetic long-range 
order. With increasing J 2 , the frustration increases, 
but at J 2 = Ji, the model again has the ground 
state with the collinear antiferromagnetic long-range or¬ 
der. The strongest frustration occurs around J 2 /J 1 = 
0.5, where nonmagnetic states such as a valence 
bond state"*’14-16,18, 22 ,29) ^ spin-liquid 

state*^’^^“^®^ have been suggested to appear, the region 
of which has recently been studied further in detail.^*’ 
The ti-t 2 -U Hubbard model with the nearest-neighbor 
(ti) and next-nearest-neighbor (^ 2 ) hopping parameters 
and the on-site repulsive interaction U has also been 
studied, where it has been shown that the critical in¬ 
teraction strength Uc of the metal-insulator transition 
increases monotonically with increasing and 

that the ground state has the Neel order at a small ^ 2/^1 


and a collinear order around ^2 = Then, the non¬ 

magnetic insulating state appears between these ordered 
states.^"*’ 

In the triangular-lattice cases, the anisotropic J-J' 
triangular Heisenberg model has been studied. In the 
isotropic case (J = J'), the 120° spiral ordered phase 
is known to be stable.In the anisotropic case, the 
Neel order is realized when J' j J is small and the spi¬ 
ral order is realized around J' j J = 1,38-49) between 
these phases, a dimer ordered phase^®^ or a spin-liquid 
phase^^’*®^ has been predicted to appear. The anisotropic 
t-t'-U triangular Hubbard model has also been stud¬ 
ied,where it has been shown that Uc increases 
with increasing t'/t: at a small t'/t a metal-insulator 
transition occurs from the metallic phase to the Neel 
ordered phase, whereas at t'/t ~ 1 a nonmagnetic in¬ 
sulating phase appears between the metallic and spiral 
ordered phases.®®’®®^ Recently, the magnetic orders in 
the triangular-lattice Heisenberg model with Ji and J 2 
have also been studied, where a nonmagnetic insulating 
phase is shown to appear between the spiral and collinear 
phases. 

In this paper, motivated by the above developments in 
the field, we study the frustrated square-lattice Hubbard 
model at half filling with the isotropic nearest-neighbor 
and anisotropic next-nearest-neighbor hopping parame¬ 
ters and clarify the metal-insulator transition, the ap¬ 
pearance of possible magnetic orderings, and the emer¬ 
gence of a nonmagnetic insulating phase. The search 
is made in a wide parameter space including square, 
crossed-square, and triangular lattices, as well as in weak 
to strong electron correlation regimes. We use the varia¬ 
tional cluster approximation (VGA) based on self-energy 
functional theory (SFT),®^“®^^ which enables us to take 
into account the quantum fluctuations of the system, so 
that we can study the effect of geometrical frustration 
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on the spin degrees of freedom and determine the crit¬ 
ical interaction strength for the spontaneous symmetry 
breaking of the model. We examine the entire regime of 
the strength of electron correlations at zero temperature, 
of which little detail is known. In particular, we compare 
our results in the strong correlation regime with those of 
the Heisenberg model, for which many studies have been 
carried out. We also compare our results with those in 
the weak correlation limit via the generalized magnetic 
susceptibility calculation and with those of the classical 
Heisenberg model calculation where the quantum spin 
fluctuations are absent. 

We will thereby show that magnetic phases with Neel, 
collinear, and spiral orders appear in relevant regions of 
the parameter space of our model and that a nonmag¬ 
netic insulating phase, caused by the quantum fluctua¬ 
tions in the frustrated spin degrees of freedom, emerges 
in a wide parameter region between the ordered phases 
obtained. The orders of the phase transitions will also be 
determined. We will summarize our results as a ground- 
state phase diagram in a full two-dimensional parameter 
space. This phase diagram will make the characterization 
of the nonmagnetic insulating phase more approachable, 
although this is beyond the scope of the present paper. 



Fig. 1. (Color online) (a) Schematic representation of the 
square-lattice Hubbard model with the isotropic nearest-neighbor 
hopping parameter fi and anisotropic next-nearest-neighbor pa¬ 
rameters t 2 and fj- (t>) Isotropic square lattice at t2 = (‘^) 

Crossed square lattice at t2 = • (<i) Isotropic triangular lat¬ 

tice at ti = t 2 and = 0- The arrows represent the directions of 
the electron spins. The sublattices are indicated by different colors. 


2. Model and Method 

We consider the frustrated Hubbard model defined on 
the two-dimensional square lattice at half filling as illus¬ 


trated in Fig. 1. The Hamiltonian is given by 

H =-t2 

(*j) {{ij)) 

((ij))' 

where cj^ is the creation operator of an electron with 
spin cr at site i and rii^ = cl^Ci^r- {i,j) indicates the 
nearest-neighbor bonds with an isotropic hopping pa¬ 
rameter ti, and {{i,j)) and {{i,j))' indicate the next- 
nearest-neighbor bonds with anisotropic hopping param¬ 
eters t 2 and t' 2 , respectively [see Fig. 1(a)]. U is the on¬ 
site Coulomb repulsion between electrons and ^ is the 
chemical potential maintaining the system at half filling. 
In the large-!/ limit, the model can be mapped onto the 
frustrated spin-1/2 Heisenberg model 

h = j,Y,s.-s, + J2 Y: s,-s, + j'2 Y 

( 2 ) 

in the second-order perturbation of the hopping parame¬ 
ters with Si = 4 a'^a/ 3 Cj,g/ 2 , where aajs is the vec¬ 

tor of Pauli matrices. The exchange coupling constants 
are given by Ji = J 2 = 4t|/C/, and J 2 = /U 

for the lattice shown in Fig. 1(a). We will compare our 
results of the Hubbard model in the strong correlation 
regime with those of the frustrated Heisenberg model, 
for which related studies have been carried out. 

We treat a wide parameter space ofO < 12 /ti < 1 
and 0 < — 1) including three limiting cases: (i) 

at t 2 = t '2 = Q [square lattice, see Fig. 1(b)], where the 
Neel order is realized, (ii) at t 2 = t '2 = ti [crossed square 
lattice, see Fig. 1(c)], where the collinear order is realized, 
and (iii) at ^2 = ti and t '2 = 0 [triangular lattice, see 
Fig. 1(d)], where the 120° spiral order is realized. We 
will calculate how the above three ordered phases change 
when the hopping parameters are varied in the ranges 
0 < t2 <ti and t) <t'2 <ti. 



Fig. 2. (Color online) Left: twelve-site square-lattice cluster used 
as a reference system in our analysis. Right: equivalent triangular- 
lattice cluster, where the three sites 1, 2, and 3 form an equilateral 
triangle. The anisotropic triangular lattice is defined as = 0 and 
tl t2- 
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We employ the VGA, which is a quantum cluster 
method based on where the grand potential 
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f2 of the original system is given by a functional of the 
self-energy. By restricting the trial self-energy to that of 
the reference system E', we obtain the grand potential 
in the thermodynamic limit as 

f7[S'] = n' + Trln(G(7^ - E')”^ - TrlnG', (3) 

where and G' are the exact grand potential and Green 
function of the reference system, respectively, and Gq 
is the noninteracting Green function. The short-range 
electron correlations within the cluster of the reference 
system are taken into account exactly. 

The advantage of the VGA is that the spontaneous 
symmetry breaking can be treated within the frame¬ 
work of the theory. Here, we introduce the Weiss fields 
for magnetic orderings as variational parameters. The 
Hamiltonian of the reference system is then given by 
H' = H + + He + Hs with 


i 

(4) 

He = h'eY,E^^--'St 

i 

(5) 

■ S^, 

(6) 


I 


where /i^, /iq, and /ig are the strengths of the Weiss fields 
for the Neel, collinear, and spiral orders, respectively. 
The wave vectors are defined as Qn = tt) for the Neel 
order and Qc = (tt, 0) or (0, tt) for the collinear order. 
For the spiral order, the unit vectors are rotated by 
120° to each other, where Oi (= 1,2,3) is the sublattice 
index of site i. The variational parameter is optimized on 
the basis of the variational principle dVl/dh' = 0 for each 
magnetic order. The solution with h' ^ 0 corresponds to 
the ordered state. 

We use the twelve-site cluster shown in Fig. 2 as the 
reference system. This cluster is convenient because we 
can treat the two-sublattice states (Neel and collinear 
states) with an equal number of up and down spins and, 
at the same time, the three-sublattice state (spiral state) 
with an equal number of three sublattice sites. Note that 
longer-period phases such as a spiral phase mentioned in 
a different system®^) and incommensurate ordered phases 
are difficult to treat in the present approach. 

3. Results of Calculations 

3.1 Strong correlation regime 

First, let us discuss the phase diagram of our model 
in the strong correlation regime U/ti = 60. The result 
is shown in Fig. 3, where the result for our Hubbard 
model in the (^ 2 /^ 1 ; ^ 2 /^ 1 ) plane as well as the same 
result converted to the Heisenberg model parameters 
(J 2 /J 1 , JU-Ji) are shown. We find three ordered phases: 
the Neel ordered phase around (^ 2 /^ 1 ,^ 2 /^ 1 ) = (OjO); ttie 
collinear ordered phase around (^ 2 /^ 1 , ^ 2 /^) = (1,1), and 
the spiral ordered phase around — (1,0) 




Fig. 3. (Color online) (a) Calculated ground-state phase dia¬ 
gram of our model at Ujtx = 60 in the (t2/ti, pls-ne and 

(b) converted phase diagram in the {J 2 IJ 2 IJi) plane. The un¬ 
colored region corresponds to the nonmagnetic insulating phase. 
The transition to the collinear phase is of the first order (or dis¬ 
continuous) and the transitions to the Neel and spiral phases are 
of the second order (or continuous). The phases along dashed lines 
(i), (ii), and (iii) shown in (a) are circumstantiated in Fig. 4. 


and (0,1). The nonmagnetic insulating phase, which 
is absent in the classical Heisenberg model (see Ap¬ 
pendix A), appears in an intermediate region between 
the three ordered phases. Thus, the quantum fluctuations 
in the frustrated spin degrees of freedom are essential 
in the emergence of the nonmagnetic insulating phase. 
As shown below, the phase transition to the collinear 
phase is of the first order (or discontinuous) and the 
phase transitions to the Neel and spiral phases are of 
the second order (or continuous). This phase diagram 
is determined on the basis of the calculated ground-state 
energies E = fl+g, (per site) and magnetic order parame¬ 
ters M (per site) defined as Mn = (2/i) J2i 
for the Neel order. Me = (2/T) (S'?) for the 

collinear order, and Ms = (2/L) ^ ■ Go. • {Si) for the 
spiral order, where ( ) stands for the ground-state expec¬ 
tation value and L is the number of sites in the system. 
In the following, we will circumstantiate the obtained 
phases, particularly along lines (i), (ii), and (iii) drawn 
in Fig. 3(a), whereby we will discuss some details of our 
calculated results in comparison with other studies. 

Along line (i): The results are shown in the left pan¬ 
els of Fig. 4, where we assume ^2 = t^. At <2 = 0, the 
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Fig. 4. (Color online) Calculated ground-state energies (upper panels) and order parameters (lower panels) for the Neel, collinear, 
spiral, and nonmagnetic insulating phases as a function of t2/4i or The left, middle, and right panels correspond to lines (i), (ii), 

and (iii) in Fig. 3(a), where we assume t2 = *2 “ t 2 = ti, respectively. The inset in (c) and (e) displays the energy difference 

between the spiral and nonmagnetic insulating phases AE, and other insets enlarge the region near the phase boundary. 


ground state is the Neel order, and with increasing t2, 
the energy of the Neel order gradually approaches the 
energy of the nonmagnetic state. At ^ 2/^1 = 0.73, the 
energy of the Neel order continuously reaches the energy 
of the nonmagnetic state and the Neel order disappears. 
The calculated order parameter indicates a continuous 
phase transition. At t 2 /ti = 1, on the other hand, the 
ground state is the collinear order. The ground-state en¬ 
ergy of the collinear order increases with decreasing t 2 , 
and at t 2 /ti = 0.79, it crosses to the nonmagnetic state, 
resulting in a discontinuous phase transition, as indicated 
by the calculated order parameter. The nonmagnetic in¬ 
sulating state thus appears at 0.73 < ^ 2/^1 < 0.79, 
which corresponds to the region 0.53 < J 2 /J 1 < 0.63 
in the Heisenberg model parameters. In comparison with 
previous studies on the J1-J2 square-lattice Heisenberg 
model, which have estimated the transition point be¬ 
tween the Neel and nonmagnetic phases to be at J2/J1 = 
0.40 — 0 . 44 ,^^’^^’32) result slightly overestimates 
the stability of the Neel order. This overestimation may 
be caused by the cluster geometry used in our calcu¬ 
lations; if we use the 2x2 site cluster as the refer¬ 
ence system, the transition occurs at J2/J1 = 0 . 42 ,^®) 
which is in good agreement with the previous studies. 
The transition point between the collinear and nonmag¬ 
netic phases, on the other hand, has been estimated to 
be at J2/J1 = 0.59 — 0.62,^’^’^^’^'^’^^) which is in good 
agreement with our result. 

Along line (ii): The results are shown in the mid¬ 


dle panels of Fig. 4, where we assume <2 = 0. With 
increasing t 2 from t 2 = 0, at which the ground state 
is the Neel order, the energy of the Neel order grad¬ 
ually approaches the energy of the nonmagnetic state, 
and at 12 /ti = 0.88, the Neel order disappears con¬ 
tinuously. The calculated order parameter indicates the 
continuous phase transition. At ^ 2 /^ = 1, on the other 
hand, the ground state is the spiral order, although the 
energy difference between the spiral and nonmagnetic 
states is very small [see the inset of Fig. 4(c)] due to 
the strong geometrical frustration of the triangular lat¬ 
tice. With decreasing ^2 from t 2 /H = 1, the ground-state 
energy of the spiral order increases gradually and ap¬ 
proaches the energy of the nonmagnetic state, and at 
t 2 /ti = 0.89, the spiral order disappears continuously, 
in agreement with the calculated order parameter. Thus, 
the nonmagnetic phase appears in a very narrow region 
of 0.88 < t 2 /ti < 0.89. The corresponding Heisenberg 
model parameters at which the Neel and spiral orders dis¬ 
appear are around J 2 I J\ = 0.79. The previous studies for 
the anisotropic triangular-lattice Heisenberg model'^^A4) 
have given values around J 2 /J 1 = 0.80 — 0.87 for the 
transition point, which are in good agreement with our 
result. 

Along line (iii): The results are shown in the right pan¬ 
els of Fig. 4, where we assume ^2 = At <2 = 0, the 
ground state is the spiral order, although the energy dif¬ 
ference from the nonmagnetic state is very small [see the 
inset of Fig. 4(e)]. With increasing ^ 2 ; the energy of the 
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Fig. 5. (Color online) Calculated ground-state phase diagrams of our model in the (£ 2 /^ 1 , ^ 2 /^ 1 ) (a) [//£i = 10, (b) f//£i = 6, 

and (c) U/ti =2. The uncolored regions in (a) and (b) correspond to the nonmagnetic insulating phase. The phases along dashed lines 
(i), (ii), and (hi) shown in (a) are circumstantiated in Fig. 6. 
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Fig. 6. (Color online) Calculated results for the order parameters M (pink, green, and red dots) of the Neel, collinear, and spiral 
phases and the single-particle gap A/£i (black crosses) at U/ti = 10 (upper panels), U/ti = 6 (middle panels), and U/ti = 2 (lower 
panels). The left, middle, and right panels correspond to the lines (i), (ii), and (hi) defined in Fig. 5(a), where we assume £2 = £ 2 ? ^'2 ~ 
and £2 = £ 1 , respectively. 


spiral order gradually approaches the energy of the non¬ 
magnetic state, and at t^/tx = 0.34, the spiral order dis¬ 
appears continuously, in agreement with the calculated 
order parameter. On the other hand, with decreasing t '2 


from t' 2 /ti = 1, at which the collinear order is stable, 
the ground-state energy of the collinear order increases 
and crosses to the nonmagnetic state at t' 2 /ti = 0.59. 
The transition is thus discontinuous, in agreement with 
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the calculated order parameter. The nonmagnetic state 
therefore appears at 0.34 < t' 2 /ti < 0.59, which cor¬ 
responds to the region 0.11 < J 2 /J 1 < 0.34 if we use 
the Heisenberg model parameters. To our knowledge, no 
comparable calculations have been made for the frus¬ 
trated Heisenberg model in this parameter region. 

3.2 Intermediate to weak correlation regime 

Next, let us discuss the phase diagram of our model in 
the intermediate to weak correlation regime. The results 
at U/ti = 10, 6, and 2 are shown in Fig. 5. The detailed 
results for the calculated single-particle gap and order 
parameters are also shown in Fig. 6 along lines (i), (ii), 
and (iii) defined above. 

At U/ti = 10, we find that the results are qualita¬ 
tively similar to those at U/ti = 60, except for the tran¬ 
sition between the Neel and spiral orders: the nonmag¬ 
netic insulating phase appears between these orders at 
U/ti = 60 but a direct first-order transition occurs at 
U/ti = 10^^’^®) with a double minimum structure in the 
grand potential. The nonmagnetic insulating phase ap¬ 
pears at 0.71 < tijtx < 0.84 along line (i), which is in 
good agreement with the previous studies,®^’where 
the values tijtx = 0.70 — 0.77 for the transition between 
the Neel and nonmagnetic phases and ti/t\ = 0.82 — 0.85 
for the transition between the collinear and nonmagnetic 
phases were reported. The nonmagnetic phase also ap¬ 
pears at 0.42 < t^jtx < 0.64 along line (iii). 

At Ujtx = 6, we find that the spiral phase disap¬ 
pears and a paramagnetic metallic phase appears in 
the triangular lattice geometry at (< 2 /^ 1 , ^ 2 /^ 1 ) — (1>0) 
or (0,1). The nonmagnetic insulating phase appears at 
0.66 < tijtx < 0.86 along line (i), the region of which 
becomes wider with decreasing value of Ujtx from 60 to 
10 and 6, which is again in good agreement with the 
previous studies.®®’The region of the nonmagnetic in¬ 
sulating phase also becomes wider along line (iii), which 
occurs between the collinear and paramagnetic metallic 
phases. Along line (ii), the nonmagnetic insulating phase 
with a small charge gap appears again, which is between 
the Neel and paramagnetic metallic phases.®®’®®^ 

At Ujtx = 2, the paramagnetic metallic phase over¬ 
whelms the collinear and nonmagnetic insulating phases, 
retaining only the Neel ordered phase around tijtx = 
t'^jtx = 0. Within the Neel phase, the charge gap opens 
only at 0 < ^ 2/^1 < 0.16 and the metallic Neel ordered 
phase appears at 0.16 < t^jtx < 0.41 along line (i) and 
at 0.16 < t^jtx < 0.47 along line (ii). The perfect Fermi 
surface nesting at t^ifix = t'^j^i = 0 and its deforma¬ 
tion away from tijtx = t'^jtx = 0 are responsible for 
these results.®^^ The generalized magnetic susceptibility 
Xo{q) calculated in the noninteracting limit of our model 
[Eq. (1)] explains this result (see Appendix B). The tran¬ 
sition between the Neel ordered metallic phase and the 
paramagnetic metallic phase is continuous along line (i) 
and discontinuous along line (ii). 


4. Summary 

We have used the VGA based on SFT to study the 
two-dimensional frustrated Hubbard model at half filling 
with the isotropic nearest-neighbor and anisotropic next- 
nearest-neighbor hopping parameters. We have particu¬ 
larly focused on the effect of geometrical frustration on 
the spin degrees of freedom of the model in a wide param¬ 
eter space including square, crossed-square, and triangu¬ 
lar lattices in a wide range of the interaction strength 
at zero temperature. We have thereby investigated the 
metal-insulator transition, the magnetic orders, and the 
emergence of the nonmagnetic insulating phase, although 
the phases with incommensurate orders or with longer- 
period orders than the cluster size used have not been 
taken into account owing to the limitation of the VGA. 
We have also calculated the ground-state phase diagram 
of the corresponding classical Heisenberg model as well 
as the generalized magnetic susceptibility in the nonin¬ 
teracting limit. 

We have thus determined the ground-state phase dia¬ 
gram of the model and found that, in the strong corre¬ 
lation regime, magnetic phases with the Neel, collinear, 
and spiral orders appear in the parameter space, and 
a nonmagnetic insulating phase, caused by the effect of 
quantum fluctuations in the frustrated spin degrees of 
freedom, emerges in the wide parameter region between 
these three ordered phases. We have also found that the 
phase transition from the Neel and spiral orders to the 
nonmagnetic phase is continuous (or a second-order tran¬ 
sition), whereas the transition from the collinear order 
to the nonmagnetic phase is discontinuous (or a first- 
order transition). We have compared our results with 
the results of the corresponding Heisenberg model calcu¬ 
lations that have been made so far and found that the 
agreement is good whenever the comparison is possible. 
We have also found that, in the intermediate correla¬ 
tion regime, the paramagnetic metallic phase begins to 
appear in the triangular lattice geometry, which over¬ 
whelms the collinear and nonmagnetic insulating phases 
in the weak correlation regime, retaining only the Neel 
ordered phase in the square lattice geometry. We hope 
that our results for the phase diagram obtained in the 
wide parameter space will encourage future studies on 
the characterization of the nonmagnetic insulating phase 
as well as on its experimental relevance. 

We thank S. Miyakoshi for useful discussions. T. K. 
acknowledges support from a JSPS Research Fellowship 
for Young Scientists. This work was supported in part 
by a Grant-in-Aid for Scientific Research (No. 26400349) 
from JSPS of Japan. 

Appendix A: Ground-State Phase Diagram of 
the Classical Heisenberg Model 

Here, we present the ground-state phase diagram of 
the classical Heisenberg model, which is defined as in 
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Fig. A*l. (Color online) Calculated ground-state phase diagram 
of the corresponding classical Heisenberg model. 


Eq. (2) but its quantum spins Si are replaced by the 
classical vectors S, so that quantum fluctuations of the 
system are completely suppressed although the frustra- 
tive features in the spin degrees of freedom are present. 
The Hamiltonian is given hy H — J{q) ' ^g in 
momentum space, where 

J{q) = Ji(cosq3, +cosqy) + J2COs{qx + qy) 

+ J^cosiq^-qy). (AT) 

The ground states of the system are calculated®®^ and 
the phase diagram is obtained as shown in Fig. AT. 
We find that the magnetically ordered ground states ap¬ 
pear in the entire parameter space examined, which in¬ 
clude the Neel order [q = ( 7 r, 7 r)], collinear order [<7 = 
(tt, 0 )], and spiral orders [q = ((?,(?) and {q',—q') with 
q = cos“^(—J 1 / 2 J 2 ) and q' = cos“^(—Ji/2J 2 )]. There¬ 
fore, comparing with the results given in the main text, 
we may conclude that the quantum fluctuations in the 
geometrically frustrated spin degrees of freedom are es¬ 
sential in the emergence of the nonmagnetic insulating 
phase discussed in the main text. 


Appendix B: Generalized Susceptibility in the 
Noninteracting Limit 

Here, we present the generalized magnetic susceptibil¬ 
ity (or Lindhard function) at zero frequency,®®^®®^ 


xo{q) = jYl 


/(cfe) - /(Cfc+q) 


(BT) 


calculated for our model [Eq. (1)] in the noninteract¬ 
ing limit, where Cfc is the corresponding noninteract¬ 
ing band disperson and /(e) is the Fermi function. The 
calculated results at temperature O.Olti are shown in 
Fig. BT, where we find that a diverging behavior appears 
only at q = ( 7 r, 7 r) in Fig. BT(a) due to perfect Fermi 
surface nesting, which yields the Neel ordered state at 
^ 2/^1 = ^ 2/^1 = 0.0 in the presence of a small but finite 



Fig. B'l. (Color online) Calculated generalized magnetic sus¬ 
ceptibility defined in Eq. (B-1) at (a) t 2 /fi = = 0-0, (b) 

* 2/^1 = 1.0, = 0.0, (c) * 2/^1 = 1.0, = 0.8, and (d) 

12 /ti = t' 2 /ti = 1.0. The corresponding Fermi surface is shown in 
each panel. 


interaction strength U. There are characteristic features 
of Xo{q) but no other diverging behaviors are found, in¬ 
dicating the absence of other magnetic orderings in the 
weak correlation limit. 
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